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Abstract — This paper presents a scheme to identify the blur 
parameters using support vector machine (SVM) Multiclass 
approach has been used to classify the length of motion blur 
and sigma parameter of atmospheric blur. Different models 
of SVM have been constructed to classify the parameters. 
Experimental results show the robustness of the proposed 
approach to classify blur parameters. 

Index Terms — Classification, support vector machine (SVM), 
motion blur, point spread function(PSF) 

I. Introduction 

The goal of image restoration is to obtain visual 
information from a degraded observation. Image restoration 
finds application in many areas like photographic blur, 
astronomical imaging, remote sensing, medical imaging etc. 
Images may be blurred by the imaging device, or the medium 
through which the light propagates. The relative motion 
between the imaging device and the object also leads to image 
blurring. In most of the applications, the blurring process is 
assumed to be linear and it is mathematically represented as 

2D convolution between the original image f(x, y) and the 
degradation function h(x,y) also known as point spread 
function (PSF). Thus g(x, y) can be written as 

= 2 -^Ci ? = - p-- y - y-- « eZ C 1 ) 
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Where * denotes the linear convolution operator and £ is 
the set of integers. Image restoration techniques are oriented 
towards approximating the true image from the degraded 
observation. In classical restoration, it is assumed that the 
PSF is known apriori and the task is just to invert the 
degradation process. However, in many practical situations, 
the point spread function [ 1] is not available and it is difficult 
to obtain the sense of the original image. This situation arises 
because of some practical constraints like characterizing the 
degradation function in astronomical imaging and may be 
due to hazard of using stronger beam for a good image 
resolution in medical imaging. In such situations blind image 
decon volution is inevitable for image restoration. There exists 
a variety of blind image restoration schemes proposed in the 
literature [8, 9, 10]. Bind image decon volution problem has 
been classified into two categories. In the first category, blur 
and its parameters are identified and the PSF is constructed. 
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Then image is restored using any of the classical 
restoration methods. In the second category PSF estimation 
and image restoration are achieved simultaneously. Our work 
falls in the first category where PSF parameters are identified 
before image restoration. In this paper, we have identified the 
length of motion blur and sigma parameter of atmospheric 
turbulence blur which decides the severity of blur. Motion 
blur parameter identification have been proposed in many 
literature[2,3,4,5,7 ]. Hough transform has been used to 
estimate blur length by detecting the lines in the spectrum of 
the degraded image [5]. Authors in [6] have proposed an 
adaptive adaline network to estimate motion length of a 
degraded image. Igor Aizenberg and his team [4] developed a 
method to restore horizontal motion blurred image. Their work 
needs a rough estimation of motion length which is achieved 
using power spectrum. Similarly in [9] a method has been 
proposed to restore images degraded with atmospheric blur. 
Principal component analysis (PCA) has been used in [8 ] to 
deal with the turbulence blur. The present work attempts to 
estimate length parameter of motion blurred image and sigma 
parameter of turbulence blurred image. We utilize the blur 
parameter subsequently to construct the PSF for restoration. 
Though there are large volume of work reported, researchers 
are still active in this direction in order to improve the image 
quality and peak signal to noise ratio (PSNR). In this paper we 
have utilized the support vector machine as a classifier to 
classify the blur parameter. 

The rest of the paper is organized as follows. Different 
blur models are described in section II. Support vector machine 
and multiclass approach has been provided in section III. 
Experimental results for the proposed approach have been 
provided in section IV. Finally section V gives the concluding 
remarks. 

II. BLUR MODELS 

A. MotionBlur Model 

The point spread function for motion blur can be described 
by the following equation. 



[O Otherwise 

where L is the length of motion blur and (j) is the angle of 
motion blur. 
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B. Atmospheric Turbulence Model 

The point spread function for atmospheric blur is described 



h(x. y )= — _ exp! 

2SKT 
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where (j controls the severity of degradation. 

HI. SVM BASED BLUR IDENTIFICATION 

A support vector machine classifies data into two classes 
by separating them into two categories. SVM is also used to 
perform multi class classification [11, 15]. It distinguishes 
one class from others by constructing a hyper plane with 
maximum margin. An optimal hyper plane is constructed by 
solving a quadratic optimization problem. SVM models are 
equivalent to the classical multi layer perceptron(MLP) neural 
network where weights are optimized in order to classify the 
test data accurately. Traditional neural networks solve a non 
convex unconstrained minimization problem. SVM models 
also use a training data set (X , Y ) like MLP and update the 
weights. However, SVM uses a kernel function and weights 
of the networks are found by solving a quadratic programming 
problem with linear constraints. SVM solves the following 
optimization problem 



min —W T W+CY{Z. 

ws>,£,f 2 



r) 



subject to y i — {w r <fi(x t ) + h) < e + 4i 
{w T <jix^ b)- y i < 
4, 4*> 0;i = l,2 5 ...jj (4) 

Where E )i is the upper training error and is the lower 

training error subject to | y, -(w r ^(x ; ) + £>)| < £ and £ is a 
threshold. C is the penalty parameter. 

A. Feature Selection 

The variance parameter decreases with the blurring 
strength in a blurred image Change of variance with blur 
length for different images are shown in Figure. 1 We have 
used variance as one of the criteria to select the samples from 
the blurred image. It is also well known that edges in an 
image play a major role for image analysis. The high frequency 
regions in a cameraman image for different motion blur length 
are shown in Figure. 2. It is evident from Figure 2 as the blur 
length increases, edges becomes smoother. The smoothness 
in the blurred image can be measured using variance. We 
construct the training set by taking the regions in the blurred 
image having higher variance. Such regions are searched by 
taking a 7x7 window and sliding this from top to bottom. If g 
represents the blurred image with blur length L, then we select 
the gray values of the window g whose variance is larger. 
The feature vector is constructed as [g (:); L] T where g (:) 
represents lexicographic ordering of the gray values of the 
window having larger variance and L is the corresponding 
blur length. 
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Figurel. Variance curve for different blurred images with blur length 

B. Blur paramter classification 

The multi class SVM has been used to classify the blur 
parameter. Several approaches to multiclass SVM have been 
proposed in the literature [12, 13, 14]. These methods 
decompose the multiclass problem into number of two class 
problems. They use different decomposition methods like 
one versus rest, one versus one etc. Using the SVM for blur 
parameter classification is straightforward. 

First we construct a training set using the feature vectors 
as described above and train the feature samples with a 
multiclass SVM. In the testing phase, features of the blurred 
image with different blur length are collected and fed to the 
trained SVM model. The output of the SVM is the identified 
blur length. The sigma parameter of the atmospheric 
turbulence blur has been identified in the same way by 
constructing another multi class SVM model. Details of the 
experiments conducted are described in the simulation 
section. 

IB. SIMULATION RESULTS AND DISCUSSION 

Experiments have been conducted in the MATLAB 
environment to show that the trained SVM model can be 
used to identify the blur parameter in new images. Standard 
images like Cameraman, Lena and Pepper are used for both 
training and testing. All the images used for training and 
testing are scaled in the range [0, 1] in the spatial domain. We 
have used a window size of 7x7. So the size of one sample 
vector used in the SVM is [49x1]. The SVM model trained 
with Cameraman image is tested with a different image. 

A. Experiment 1 : Cameraman SVM model for motion blurred 
image 

The Cameraman image is horizontally motion blurred with 
L = 2 to 30 with a gap of 2. The 1 844 samples from the blurred 
image have been taken by considering the window whose 
variance is larger than 0.02. The SVM istrained with the 
collected samples of blurred cameraman image and a SVM 
model is created. 
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Figure 2 Behavior of an edge under motion blur of various Blur length with 45° angle 



The SVM model is validated by testing the samples of 
cameraman image as well as Lena image. The vote count for 
each of the estimated class is for cameraman image and Lena 
image are noted and has been shown in TABLES I & II 
respectively. Values in each column show the number of 
votes won by different blur lengths. The maximum vote 
obtained by any row in a particular column is the estimated 
class which is equal to the blur length We should note that 
the samples from Lena image are different from the sample 
used in the training. It is clear from the Table 1 that the 
proposed multiclass SVM identifies the blur length accurately 
for the same image as well as different image. To create the 
SVM model LIBSVM package [16] is used with the default 
radial basis function. The value of C is chosen on a trial and 
error basis and for the cameraman image SVM model it is 
taken as 1 . Error tolerance and Gamma value are chosen to be 
0.0001 and 5 respectively. 

TABLEI. Blur Idetification Performanceof Multiclass Svm For 
Cameraman Image 



Actual Blur Length 



1 





913 


480 


276 


175 




5 


9 


13 


17 


5 


S61 


128 


29 


16 


9 


44 


321 


64 


28 


13 


3 


22 


154 


41 


17 


5 


9 


30 


88 



TABLE II. Blur Idetification Performanceof Multiclass Svm For Lena 
Image 



Actual Blur Length 



Id 





756 
5 


515 
9 


326 
13 


190 
17 


5 


689 


156 


51 


19 


9 


53 


331 


72 


35 


13 


6 


16 


168 


39 


17 


8 


12 


35 


97 



These parameters can also be optimized for different 
images and different degradation. After the blur parameter 
is identified PSF is constructed and the degrade images are 
restored. The restoration result for the cameraman image is 
shown in Figure 3. 



B. Experiment 2 Pepper SVM model for atmospheric blurred 

IMAGE 

The Pepper image is degraded with atmospheric turbulence 
with different sigma (a) values in the range [0.5 -5] using the 
inbuilt function available in MATLAB. The training samples 
of the blurred Pepper image are accumulated using the same 
procedure as explained in Experiment 1 Total 1392 samples; 
have been used for training the SVM.. To validate the 
effectiveness of the proposed scheme, testing samples from 
cameraman image also have been used. The performance of 
the multi class SVM for different images has been shown in 
TABLE III and IV respectively. The results in the TABLE'S 
show that the proposed multiclass SVM also works well for 
images degraded with atmospheric turbulence blur and it 
achieves generalization in terms of different images. Same 
voting strategy has been followed in this experiment to 
identify the actual sigma value. The value of C has been kept 
1 and RBF kernel has been used in the LIBSVM package. 
Gamma value and error tolerance are chosen to be 4 and 0.0001 
respectively. The restored image is obtained after constructing 
the PSF from the identified sigma parameter. The restoration 
results for Pepper image is shown in Figure4. 

TABLE III. Blur Idetification Performanceof Multiclass Svm For 
Pepper Image 
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TABLE IV. BLUR IDETIFICATION PERFORMANCEOF MULTICLASS SVM FOR 
CAMERAMAN IMAGE 
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CONCLUSIONS 

A novel method has been proposed to identify the blur 
parameter from the blurred images. The blur identification 
problem is modeled as a multiclass classification problem 
and has been solved using support vector machine. The 
suggested scheme achieves generalization in terms of 
different images and estimates the blur parameters accurately. 
The proposed scheme is applicable to images degraded with 
horizontal motion blur and images degraded with atmospheric 
turbulence blur. However, the proposed method can also be 
applied to estimate angle of motion blur by constructing 
another SVM model. 
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Figure 3. (a) Original Image (b) B Quired image with L = 17 (c ) Restored image after 
parameter identification 
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Figure 4 _{a) Original Image 


(b) Blurred image with sigma 


= 3 


(c ) Restored image after 
parameter identification 
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